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Abstract 

A simple and analytical formula is suggested to solve the problems of the 
local burnup and the isotope distributions. The present method considers 
two extreme conditions of neutrons penetrating the fuel rod. Based on these 
considerations, the formula is obtained to calculate the reaction rates of 

and ^^®Pu and straightforward the local burnup and the isotope distri¬ 
butions. Starting from an initial burnup level, the parameters of the formula 
are htted to the reaction rates given by a Monte Carlo (MC) calculation. 
Then the present formula independently gives very similar results as the MC 
calculation from the starting to high burnup level, but takes just a few min¬ 
utes. The relative reaction rates are found to be almost independent on the 
radius (except (n, 7) of and the burnup, providing a solid background 
for the present formula. A more realistic examination is also performed when 
the fuel rods locate in an assembly. A combination of the present formula 
and the MC calculation is expected to have a nice balance on the accuracy 
and the cost on time. 


1. Introduction 

To increase the efficiency of the fuel, one possible way is to increase the 
burnup of the fuel before discharge. The local burnup on the edge of the 
UO2 fuel rod is much higher than the average burnup. Thus it is of great 
importance to investigate the properties of the rim of the fuel rod when 
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considering the increment of the average bnrnnp. Many investigations show 
that the mechanical structnre close to the snrface is rather different from 
that at the center of the fnel rod. In the high bnrnnp range, a microstrnc- 
tnre change was fonnd on the rim of the fnel rod throngh transmission elec¬ 
tron microscopy (H- One explanation of the formation of the high bnrnnp 
strnctnre snpposed that the babbles in high bnrnnp region are nncleated and 
stabilized by hssion fragments, which depends on the hssion rate j^. The 
small pores at the high bnrnnp region are calcnlated to be highly overpres- 
snrized [^. Recently, the UO2 fnel in both light and heavy water reactor are 
investigated to nnderstand the high bnrnnp strnctnre of the fnel BB- 

In a thermal reactor, the nentrons generated from hssion need to be slowed 
down in the moderator to be able to indnce next hssion. When the low energy 
nentrons go from the moderator to the fnel rod, they are hrstly absorbed by 
the fnel close to the snrface. In general, the reaction rate is higher in the 
rim when indnced by slowed nentrons, snch as the (n, /) reaction of ^sspj 
and ^^®Pn indnced mainly by the thermal nentron. In contrast, the (n,7) 
reaction of jg j^iainly indnced by the resonance nentrons. The local 
bnrnnp phenomena is mainly cansed by the (n, 7) reaction of which 
has large cross section for resonance nentrons, especially some high peaks of 
cross section at certain energies. Snch reaction prodnces mnch more ^^®Pn 
near the snrface of the fnel rod, thns giving origin to the increased local 
bnrnnp. As a consequence of the different reaction rates along the radial 
direction, the U and Pu isotopes also have different radial distributions. In a 
fast reactor, the Pu isotopes are redistributed along the radial direction B- 
A large local bnrnnp is not expected because the neutrons do not need to be 
slowed down in the moderator. 

In nuclear science, theoretical models are very important and powerful 
to solve various of problems. Our previous works have shown that the mass, 
levels, electromagnetic moments and transitions, and other properties of the 
nuclei can be described with high precision in the frame of the nuclear shell 
model even for very unstable nuclei, such as and ^^Al B- nuclear reac¬ 
tor, the simulation is even more important because the experimental data are 
relatively difficult to be obtained, sometimes with high risk. For the present 
problem, various models can be used to calculate the local burn up and 
the isotope distributions along the radial direction. The TRANSURANUS 
model is very successful in the description of the local bnrnnp and related 
properties in various reactors including light and heavy water reactors js-lO 


In TRANSURANUS model, a fixed neutron flux is assumed except for the 
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absorption reaction of while the RAPID model consid 


ers detailed properties of the neutron flux at each burnup H- Recently, 
an empirical formulation is suggested to describe the local burnup and high 
burnup properties based on DIONISIO code 
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Here we suggest an analytical and simple formula to calculate local bur¬ 
nup and isotope distributions based on the constant relative reaction rates 
at the different radius (except for and burnup. The present method 

considers the radial distribution of the neutron flux through two extreme 
conditions. The parameters are htted to the reaction rates of the starting 
burnup level. From this level, the present method are shown to have inde¬ 
pendent and similar results as the Monte Carlo (MC) simulation, but just a 
few minutes is needed because its simplicity. 


2. The description of the model 

The rate of certain reaction type (x) between neutrons and a given nuclide 
(A) atan average burnup (bu) can be written as a function of the radius (r) 
and the energy of neutron (B). 

R(r,A,x) = I N{r, A)(j){r, E)a{A,x, E)dE, ( 1 ) 

JE 

where N (r. A) is the concentration of the nuclide A at the radial position 
r, (j){r,E) is the neutron flux at the energy E and the radius r, a{A,x,E) 
is the microscopic cross section between neutron and A for reaction type x 
at the energy E. The (n, /) and (n, 7) reactions are assumed to be the two 
more important reaction types in the fuel rod. In the present work, the main 
purpose is to examine the validity of the present model. Thus only the 

and isotopes are considered to constrain the model testing to a 

simple case. 

In the present study, the Eq. ([T]) is assumed to be two terms for (n, /) 
and (n, 7) of and and the (n, /) of 

R{r) = N{r)[ci + C2e^^^^-'^% ( 2 ) 

and three terms for (n,7) of 

R(r) = Ar(r)[c3 + C4e''i("-"“)-l-C5e^2(’'-'’'’)], ( 3 ) 
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where the coefficients c are specihed for each nuclide and the reaction type. 
The coefficient hi is expected to be the same for all three nuclides. The 
coefficient 62 in the (n, 7) of jg much larger than hi. 

The following discussion concentrates on the explanations of the above 
assumptions of the reaction rates. In a thermal reactor, two kinds of the neu¬ 
tron reactions are most important, the scattering reaction in the moderator 
and the absorption reaction in the fuel rod. When the slowed down neutrons 
penetrate into the fuel rod, we mainly concentrate on the neutron flux at the 
thermal and resonance region which are important to the absorption reac¬ 
tion in the fuel rod. The average energy of hssion neutron is around 2 MeV, 
which is much larger than that of the thermal and resonance neutron. Thus 
no multiplication of the neutron flux is considered at the energy of thermal 
and resonance region in the following discussion. Two extreme conditions 
can be considered for certain energy of the neutron. The hrst one is that the 
neutrons are strongly absorbed by the fuel, which indicates that the mean 
free path (MFP) of the absorption, A(a,E) = ^ = s.,5(i;)+s..8(i;)+s.,9(i;) ^ 
is much smaller than the size of the fuel rod at the energy E. The second 
one is that few neutrons are absorbed in the fuel at the energy E. 

For the hrst extreme condition, one can further simplify that the velocity 
of the neutrons are all perpendicular to the surface of the fuel rod. Such 
assumption is acceptable because the velocity is almost symmetric around the 
radial directions in the real case. Then the transportation of the neutrons at 
certain energy E is the solution of the Boltzmann equation in the cylindrical 
coordinate without the terms of the scattering and the source. The radial 
part of the equation is: 


1 ( 9 r 0 (r, E) 
r dr 


+ S„(E)0(r,E) 


0 , 


( 4 ) 


where the negative sign comes from the opposite direction between the ve¬ 
locity and the radial vector, and the solution is: 


(j){r,E) = ( 5 ) 

r 


where tq is the radius of the fuel rod. Because of the assumption of the 
strong absorption, the neutron flux is meaningful near the surface of the fuel 
rod. The divergence at r = 0 should not be considered. Because \{a,E) is 
much smaller than the size of the fuel rod, tq, the neutron flux decreases very 
quickly to zero at small (r — tq). The neutron flux E) is approximately 
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0(ro, near the surface and zero at the other radial region. The 

reaction rate per nuclei i?(r, A, x, E)/N{r, ^4) is proportional to 
resulting the third term in Eq.(|3]). One example of such situation corresponds 
to the {n, 7 ) reaction of certain energy. The magnitude of the atomic 

concentration of jg ^j-ound 10^^/cm^ in the fuel rod. For (n, 7 ) reaction 
of some peaks of cross section at the resonance region can achieve 10 ^ 
barn [l5|, resulting A(a, E) = ^(e)Ns magnitude of 10“^cm, which is 

much smaller than the size of the fuel rod. 

For the second extreme condition, the X{a,E) is much larger than the 
size of the fuel rod tq. The 0(r, E) is almost unchanged in the fuel rod 
and R{r, A,x, E)/N{r, A) is also approximately constant, resulting the hrst 
term in Fq.([2]) and (j3]). The magnitude of the concentrations of and 
are around or less than 10^°/cm^ in the fuel rod. If the (Ta^ 5 -a, 9 {E) and 
the aa^siE) are much smaller than 10 ^ barn and 10 ^ barn, respectively, the 
X{a,E) is much larger than tq. 

For the situations between these two limits, the second term in Fq.fl2l) 
and ([3]) are assumed with 61 the same magnitude of Please note that 
the above discussion is restricted in the fuel rod. If the neutrons go out of 
the fuel rod, they may be slowed and re-enter the fuel rod with a different 
energy. 

In principle, c, 6 in Fq.dl]) and ([3]) can be calculated through the cross 
section data and the neutron flux at tq. But in some energy regions, the 
cross sections changes dramatically, and hence it is difficult to obtain the 
coefficients. Spatially Dependent Dancoff Method can be used to calculated 


the cross-sections in the resonance region [16(]. In the present work, these 
coefficients are htted using a MC calculation. An initial burnup level bui is 
assumed and with it a MC calculation is done for one fuel rod to calculate 
the reaction rates at different radial positions. The coefficients are then 
htted to these reaction rates and used to calculate the burnup levels after 
bui- After bui, the MC calculation and the present formula are performed 
independently for all burnup levels, bu 2 , bu^ and so on. The concentrations 
can be calculated through: 


A^5(r, bui+i) = A^5(r, bui) - {Rnjir) -h R5^^{r))T, 

Ns{r, bui+i) = Ns{r, bui) - {Rgjir) Rs,^ir))T, 

A^9(r, bui+i) = A^9(r, bm) + {Rs,'y{r) - R9jir))T, 


( 6 ) 


by assuming that the reaction rates do not change during the time T. The 
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time duration T corresponds to the burnup change between two levels: 


A A + Rsjir) + RQj{r)]27irdr 

-—-7--. 

Mo[/(7rr2) 

where Q is the average energy released by hssion and MQu{nrQ) is the initial 
mass of U isotopes in volume TTr^. One can transform: 


AN,{r) 

^Ns{r) 

ANg{r) 


Mpujirrl) _-(^5j(r) +^5,7(^))_ 

Q /J’°[^ 5 ,/(r) + R%j{r) + Rgj{r)]2nrdr 
Mou{Tirl) _ -{R^j{r) + R^^^{r)) _ 

Q /J’°[^ 5 ,/(r) + R^j{r) + i?9j(r)]27rrdr 
Mpujirrl) _ (i?8,7(^) - R9,f{^)) _A, 

Q Io°[R5,f{r) + Rsjir) + Rgj{r)]27irdr 


In MC calculation, the reaction rates in above equations are simulated with 
the concentration of all three isotopes at each burnup level. The present for¬ 
mula calculate all reaction rates through the coefficient fitted to the reaction 
rate at bui. The Rp^^{r) reaction rate is neglected in above equations in both 
MC and analytical calculation because the present model does not include 
^^°Pu. It is acceptable as present work concentrates on the examination of 
the formula not on a real burnup problem. The next section shows that the 
present formula can obtain results for bun quite close to those obtained with 
the MC calculation. 

The above equation looks similar to the formula in TRANSURANUS 
model js, 10|, in which the relationship between N{r) and bu can be generally 
written as: 


dNj(r) 

dbu 






( 8 ) 


where j is one kind of nuclide. For 235,238u^ there are no the second terms. 
For 236U, 237Np^ and 238,239,240,241,242p^^ ^ ^ 235,236u^ 237Np^ 238u^ 

and 239,240,241 respectively. The a a and cxc are the one-group effective cross 
sections for the total neutron absorption and neutron captured, respectively. 
The cross sections are obtain differently for UO 2 and MOX fuel because of 
the very different initial concentrations in each fuel and the corresponding 
different neutron spectrum. A is a conversion constant. The /(r) is the radial 
form factor, which is 1 + piexp{—p 2 {ro — ry^) for 238p 240pg 
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Table 1: Relative radius, atomic concentrations and relative reaction rates of each nucleus. 
The atomic concentrations are given in cm“^, all the other quantities are dimensionless. 


r'/ro 

N 5 

Ns 

No 

Rkf 

Rkf 

Rkf 

Rk-y 

R's.y 

1/4 

6.65x10^*^ 

2.17x1022 

3.69x10^0 

1.00 

0.0025 

2.50 

0.23 

0.0183 

5/8 

6.63x1020 

2.17x1022 

3.91x1010 

1.02 

0.0024 

2.49 

0.23 

0.0189 

13/16 

6.61x1020 

2.17x1022 

4.33x1010 

1.05 

0.0024 

2.49 

0.23 

0.0203 

29/32 

6.59x1020 

2.17x1022 

5.01x1010 

1.07 

0.0023 

2.48 

0.23 

0.0230 

61/64 

6.58x1020 

2.17x1022 

6.26x1010 

1.08 

0.0023 

2.47 

0.23 

0.0283 

125/128 

6.58x1020 

2.16x1022 

8.49x1010 

1.09 

0.0022 

2.47 

0.23 

0.0368 

127/128 

6.57x1020 

2.16x1022 

1.26x1020 

1.10 

0.0022 

2.47 

0.23 

0.0552 


for all other nuclide. The /(r) comes from the resonance absorption and the 
parameters are determined by comparison with measurements The 

local burnup and isotope distributions can be calculated through Eq. (|8]). 
More details can be found in Ref. j^, [l^. The present work considers not 
unit /(r) (actually the different radial neutron flux) for all nuclei. 
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3. Calculations and discussions 

In the present work, the continuous energy Monte Carlo code TRIPOLI- 
is used for the MC calculations as the starting point of the analytical 
formula and the reference for after calculations. The geometry of the fuel 
rod in the present investigation is set to be tq = 0.4127 cm, with cladding 
between 0.4127 and 0.4744 cm. The moderator is in an outside box with 
the length 1.2647 cm with all surface reflection all neutrons. The fuel rod is 
divided to seven parts in the MC simulation, with the dividing point located 
at r/ro = 1/2,3/4,7/8,15/16,31/32,63/64. The corresponding center r'/ro 
of each part is listed in Table [TJ For a normal UO 2 fuel, there is no Pu 
isotopes at the beginning. It is reasonable to start present calculation at a 
certain burnup bui with Pu included. The MC calculation is done with 3.3% 
enrichment f^gj ggjj obtain [n, f) and (n, 7 ) reaction rates 

of and With the reaction rates, the number of and 

^^®Pu in a certain burnup can be calculated by assuming the reaction rates 
do not change in this period. The concentration of and ^^®Pu at 

bui = 3.4 MWd/kgU is given in Table [T]as the starting point of the following 
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calculations. 

The (n, /) and (n, 7 ) reaction rates of and can be sim¬ 

ulated by using concentrations listed in Table [H The corresponding relative 
reaction rates defined below are also listed in Table [1] 


KA’-) 

Kf(r) 

n's-A') 


R<iAr)mr) 

R^jir = ro/4)/Af5(r = ra/ij) 

RsAA/NsA) 

R^AAmA 

It,j(r)/N,(r) 

R^AA/mA 

R^iyA) 

R<i.f(A 

R,AA/na-) 

R^AA/mA 


The above equations firstly scale the reaction rate by the corresponding con¬ 
centration. Because the concentration is different at different radial position. 
The scaling cancels such difference in the consideration of reaction rate. And 
then calculate the reaction rates relative to that of {n, f) reaction of We 
hnd that the radial distribution of these relative reaction rates are almost 
constant except {n, 7 ) reaction of Thus all the reaction rates except 

that of (n, 7 ) of pg treated in Eq ([2]) with the same coefficients. 

The radial distributions of these reaction rates are fitted as, 


R 5 j{r) = cW(r)/(r), 

Rgjir) = 0.024:cN8{r)f{r), 

Rgjir) = 2A9cNg{r)f{r), 

R5,'yir) = 0.23cN5{r)f{r), 

/(r) = + 

/?8,7(r) = 0.0183cA'8(r’)[0.97 + 0.42e^-^(’'/'’“-^)(9) 


where c is the reaction rate per nuclei of {n, f) of ^ _ j-q/A and 

canceled in Eq. d?]). Please note that the same bi is used for all nuclei. From 
bu 2 , MC and analytical calculations are independent. It is straightforward to 
see the advantage of transforming the burnup problem to the reaction rate 
problem in the present work. The reaction rate actually is the most direct 








Radius (r/r^) 


Figure 1: The comparison of the local burnup given by present formula and MC calcula¬ 
tions. 

quantity to obtain the burnup and the concentration. In many models, the 
reaction rate must be obtained after solving the neutron flux problem. The 
present work treats the neutron flux problem in a simple way in section [2] 
and considers the relative reaction rate. The validity of these considerations 
will be discussed later. 

In the following discussion Mm is set to be 2.85 MWd/kgU and the total 
burnup is calculated to hui^ = 43.30 MWd/kgU. The local burnup from 
two calculations are presented in FIG. [1] The analytical calculations can 
give almost the same results as MC calculations. The detailed data of three 
burnup leves in Table. [2] shows that the average errors are less than 1%. 

The radial distributions of U and Pu isotopes at three burunp levels are 
presented in Table |3l The bu 2 is the hrst burnup level that MC calculation 
and present formula are performed independently. The error between two 
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Table 2: The comparison of the local burnup (unit in MWd/kgU) at three levels given by 
present formula (AN) and MC calculations. 


r'/ro 

bU2,MC 

bU2,AN 

bUs,MC 

bUs^AN 

bUu^MC 

bUu^AN 

1/4 

5.96 

5.94 

21.48 

21.33 

36.57 

36.30 

5/8 

6.11 

6.09 

22.11 

22.09 

37.69 

37.80 

13/16 

6.29 

6.29 

23.05 

23.11 

39.53 

39.75 

29/32 

6.48 

6.49 

24.40 

24.22 

42.42 

41.83 

61/64 

6.72 

6.74 

26.62 

26.37 

47.48 

46.55 

125/128 

7.06 

7.08 

30.13 

30.38 

55.50 

55.91 

127/128 

7.68 

7.72 

37.20 

37.31 

71.75 

71.64 

average errors 


0.32% 


0.55% 


0.84% 


methods at this level reflects how exactly the Eq. ([9]) fits to the MC reaction 
rates. The concentration of changes around 10% from hui (in Table [1]) 
to bu 2 (in Table E]). The average error 0.05% indicates that the present 
formula well describes the isotope distribution of at the burnup level. 
The average error of the isotope distribution of ^^®Pu is around 1%. It is 
partially from the large change of the concentration from hui to bu 2 and 
partially from the coefficient of the (n, /) reaction rate for ^^®Pu in Eq. (|9]). 
The coefficient in Table [U is decreasing from 2.50 to 2.47. The use of 2.49 
in the Eq. [9] contributes to the 1% error between MC calculation and the 
present formula at bu 2 - 

At the burnup level bu^ and 5 m 14 , the errors become larger, indicating that 
the reaction rates in Eq. ([9]) change slightly during burnup. The detail will 
be discussed later. It can be seen that the radial distribution of Ng does not 
change much after a few periods because the production and reaction rates 
find a balance at such distribution. Both the comparisons of the local burnup 
and the isotope distributions imply that the present simple and analytical 
formula can give very nice description if a starting point is given. 

As the present model is constrained in a single rod, it is worth to examine 
its validity in a more realistic case, such as in an assembly. Figure | 2 ] shows the 
model of an assembly which has the same size for each fuel rod as the previous 
single one. There are 24 positions for the control rods and one postion in the 
center for the detectors. In the present calculation, the control rods are not 
inserted and the rod for detectors are filled with water. The reaction rates of 
the fuel rods in the assembly are simulated by setting all fuel rods the same 
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Table 3: The comparison of the concentrations (unit in cm at three burnup levels given 
by present formula (AN) and MC calculations. 


r'/ro 


Ns,mc 

Ng^MC 

-^5,TV 


Ng^AN 

1/4 

bU2 

6.01E+20 

2.17E+22 

6.13E+19 

6.00E+20 

2.17E+22 

6.09E+19 

5/8 

5.98E+20 

2.17E+22 

6.46E+19 

5.98E+20 

2.16E+22 

6.43E+19 

13/16 

5.94E+20 

2.16E+22 

7.09E+19 

5.94E+20 

2.16E+22 

7.08E+19 

29/32 

5.91E+20 

2.16E+22 

8.13E+19 

5.91E+20 

2.16E+22 

8.16E+19 

61/64 

5.89E+20 

2.16E+22 

l.OlE+20 

5.89E+20 

2.16E+22 

1.02E+20 

125/128 

5.88E+20 

2.16E+22 

1.31E+20 

5.88E+20 

2.16E+22 

1.35E+20 

127/128 

5.87E+20 

2.15E+22 

1.98E+20 

5.87E+20 

2.15E+22 

2.02E+20 

average errors 




0.05% 

0.06% 

1.06% 

1/4 

bug 

3.20E+20 

2.14E+22 

1.42E+20 

3.24E+20 

2.14E+22 

1.33E+20 

5/8 

3.14E+20 

2.14E+22 

1.49E+20 

3.18E+20 

2.14E+22 

1.42E+20 

13/16 

3.07E+20 

2.14E+22 

1.61E+20 

3.10E+20 

2.14E+22 

1.53E+20 

29/32 

3.01E+20 

2.13E+22 

1.84E+20 

3.04E+20 

2.14E+22 

1.65E+20 

61/64 

2.98E+20 

2.13E+22 

2.26E+20 

3.01E+20 

2.13E+22 

2.02E+20 

125/128 

2.96E+20 

2.11E+22 

2.93E+20 

2.99E+20 

2.11E+22 

2.78E+20 

127/128 

2.94E+20 

2.08E+22 

4.36E+20 

2.97E+20 

2.09E+22 

4.06E+20 

average errors 




1.07% 

0.03% 

7.03% 

1/4 

buu 

1.53E+20 

2.12E+22 

1.65E+20 

1.60E+20 

2.12E+22 

1.52E+20 

5/8 

1.48E+20 

2.12E+22 

1.70E+20 

1.54E+20 

2.12E+22 

1.61E+20 

13/16 

1.42E+20 

2.11E+22 

1.83E+20 

1.47E+20 

2.11E+22 

1.72E+20 

29/32 

1.38E+20 

2.10E+22 

2.07E+20 

1.42E+20 

2.11E+22 

1.85E+20 

61/64 

1.35E+20 

2.09E+22 

2.54E+20 

1.39E+20 

2.09E+22 

2.23E+20 

125/128 

1.33E+20 

2.06E+22 

3.27E+20 

1.37E+20 

2.06E+22 

3.07E+20 

127/128 

1.32E+20 

2.01E+22 

4.79E+20 

1.36E+20 

2.01E+22 

4.43E+20 

average errors 




3.36% 

0.07% 

7.89% 
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Figure 2: The model of an assembly, with filled and empty circle for fuel and control rods, 
respectively. Four selected rods are marked. 
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concentrations as the single rod at a certain burnup level. The single rod 
model are compared with four selected fuel rods in the assmebly, marked in 

FIG.H 

The relative reaction rates, dehned by Eq. ([9]), in the single rod and four 
rods in the assembly are compared at the burnup level bu 2 in Table HI The 
relative reaction rates of all hve rods are generally similar, especially for 
the two important reactions, {n, f) reaction of and The relative 

reaction rate of {n, 7 ) reaction of in the single rod is a few percent less 
than those in the four rods in the assembly. A slight modihcation can be 
applied for the Eq. ([9]) to £t to the fuel rods in the assembly. In the case of 
assembly, it is clear to see that the relative reaction rates are almost constant 
for all four fuel rods, while the reaction rates are a little different depending 
on the position of the fuel rods, seen from in Table HI 

Table [5] presents the comparison of the relative reaction rates between 
the single rod and the rodl in the assembly at the burnup level bug and 
6 ^ 15 . As relative reaction rates in the four rods in the assembly are similar 
to each other, only the data of rodl are shown in Table [5l At both bug and 
bui5, the {n, f) reaction of similar in both rods, while 

the relative reaction rate of {n, 7 ) of reaction single rod is a few 

percents larger than that in the rodl. The situations are similar to that at 
bu 2 - The comparisons in Table HI and [5] imply that the single rod model is a 
nice estimation of the fuel rods in the assembly. 

It is worth to consider why such simple formula Eq. ([9]) works. Take 
R'gj{r) for example. 



a<,Ar)IN^(r) 

/e E)a(9,E)dE 
J^4>{r,E)a(5,E)dE' 


which is the ratio of the one-group effective cross section. An approximately 
constant R'g ^{r) indicate a constant neutron dynamics for the ratio between 
the hssion reaction of and 235pj along the radial direction. Although 

the one-group effective cross section changes as the function of r, the ratio of 
the one-group effective cross section of each reaction included in the present 
study keeps almost the same except the (n, 7 ) reaction of One of the 
advantage of the present method is the consideration of the reaction rate 

[iMi 


compared with the previous methods 


Because the local burnup and 
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Table 4: The comparison of the relative reaction rates for a single fuel rod and four selected 
fuel rods in an assembly at bu 2 - 


r '/tq 


7a 

00 


-^5,7 

7a 

00 

R^j{rodi) 
Rg.f (rodl) 


single 






1/4 

1.00 

0.0026 

2.43 

0.23 

0.0186 


5/8 

1.02 

0.0025 

2.43 

0.23 

0.0192 


13/16 

1.05 

0.0024 

2.43 

0.23 

0.0206 


29/32 

1.07 

0.0024 

2.42 

0.23 

0.0235 


61/64 

1.08 

0.0023 

2.42 

0.23 

0.0289 


125/128 

1.09 

0.0023 

2.42 

0.23 

0.0376 


127/128 

1.09 

0.0023 

2.42 

0.23 

0.0561 



rodl 






1/4 

1.00 

0.0025 

2.44 

0.23 

0.0178 


5/8 

1.02 

0.0024 

2.43 

0.23 

0.0185 


13/16 

1.05 

0.0023 

2.42 

0.23 

0.0197 


29/32 

1.07 

0.0023 

2.42 

0.23 

0.0222 


61/64 

1.08 

0.0022 

2.42 

0.23 

0.0269 


125/128 

1.09 

0.0022 

2.42 

0.23 

0.0352 


127/128 

1.10 

0.0022 

2.42 

0.23 

0.0521 



rod2 






1/4 

1.00 

0.0025 

2.42 

0.23 

0.0179 

0.99 

5/8 

1.02 

0.0024 

2.42 

0.23 

0.0184 

0.99 

13/16 

1.05 

0.0023 

2.41 

0.23 

0.0197 

1.00 

29/32 

1.07 

0.0023 

2.40 

0.23 

0.0225 

1.00 

61/64 

1.09 

0.0022 

2.40 

0.23 

0.0277 

1.00 

125/128 

1.09 

0.0022 

2.40 

0.23 

0.0355 

1.00 

127/128 

1.10 

0.0022 

2.40 

0.23 

0.0538 

1.00 


rod3 






1/4 

1.00 

0.0023 

2.43 

0.23 

0.0173 

1.05 

5/8 

1.02 

0.0023 

2.43 

0.23 

0.0178 

1.05 

13/16 

1.05 

0.0022 

2.42 

0.22 

0.0191 

1.06 

29/32 

1.07 

0.0021 

2.41 

0.22 

0.0220 

1.06 

61/64 

1.09 

0.0021 

2.41 

0.22 

0.0271 

1.06 

125/128 

1.09 

0.0021 

2.41 

0.22 

0.0358 

1.06 

127/128 

1.10 

0.0021 

2.41 

0.22 

0.0533 

1.06 


rod4 






1/4 

1.00 

0.0023 

2.42 

0.22 

0.0170 

1.08 

5/8 

1.03 

0.0022 

2.42 

0.22 

0.0175 

1.09 

13/16 

1.05 

0.0021 

W 

0.22 

0.0190 

1.09 

29/32 

1.07 

0.0021 

2.40 

0.22 

0.0216 

1.09 

61/64 

1.08 

0.0020 

2.40 

0.22 

0.0268 

1.09 

125/128 

1.09 

0.0020 

2.40 

0.22 

0.0349 

1.09 

127/128 

1.10 

0.0020 

2.39 

0.22 

0.0526 

1.09 



Table 5: The comparison of the relative reaction rates for a single fuel rod and the fuel 
rod in an assembly at bug and buig. 


r '/tq 

Rkf 



K, 

K. 

1/4 

bug single 
1.00 

0.0025 

2.23 

0.23 

0.0183 

5/8 

1.02 

0.0025 

2.24 

0.23 

0.0190 

13/16 

1.04 

0.0024 

2.25 

0.23 

0.0205 

29/32 

1.06 

0.0023 

2.26 

0.23 

0.0233 

61/64 

1.07 

0.0023 

2.26 

0.23 

0.0287 

125/128 

1.08 

0.0023 

2.27 

0.23 

0.0374 

127/128 

1.08 

0.0023 

2.27 

0.23 

0.0559 

1/4 

bug rodl 
1.00 

0.0024 

2.23 

0.23 

0.0175 

5/8 

1.02 

0.0024 

2.24 

0.23 

0.0182 

13/16 

1.04 

0.0023 

2.25 

0.23 

0.0197 

29/32 

1.06 

0.0022 

2.26 

0.23 

0.0222 

61/64 

1.07 

0.0022 

2.26 

0.23 

0.0272 

125/128 

1.07 

0.0022 

2.27 

0.23 

0.0354 

127/128 

1.08 

0.0022 

2.28 

0.23 

0.0523 

1/4 

bui 5 single 
1.00 

0.0022 

2.18 

0.23 

0.0168 

5/8 

1.02 

0.0022 

2.19 

0.23 

0.0174 

13/16 

1.04 

0.0021 

2.20 

0.22 

0.0188 

29/32 

1.05 

0.0021 

2.21 

0.22 

0.0214 

61/64 

1.06 

0.0021 

2.22 

0.22 

0.0265 

125/128 

1.06 

0.0020 

2.22 

0.23 

0.0346 

127/128 

1.07 

0.0020 

2.23 

0.23 

0.0510 

1/4 

6 -Ui 5 rodl 
1.00 

0.0021 

2.18 

0.22 

0.0161 

5/8 

1.02 

0.0021 

2.19 

0.22 

0.0166 

13/16 

1.04 

0.0020 

2.20 

0.22 

0.0178 

29/32 

1.05 

0.0020 

2.20 

0.22 

0.0203 

61/64 

1.06 

0.0020 

2.22 

0.22 

0.0251 

125/128 

1.06 

0.0020 

2.21 

0.22 

0.0327 

127/128 

1.07 

0.0019 

2.22 

0.22 

0.0486 


15 



the isotope distributions do not directly depend on the neutron flux but on 
the reaction rate, it is an alternative method which consider the properties of 
the latter one. The TRANSURANUS model considers the constant neutron 
flux and the one-group effective cross sections during Abu. The assumption 
in the present model is actually the constant ratio of the one-group effective 
cross sections, which is relatively easy to be achieved. 

From the Table [H 01 and 01 it is seen that the relative reaction rates 
changes not much along the radial direction and among each burnup level. 
It should be mentioned that the relative {n, 7 ) reaction rates of ^^®Pu are also 
close to constant. It is neglected in the present calculation because it links to 
other Pu isotopes which are not included in the present study. If the error of a 
few percents is acceptable, the Eq. (0]) can be performed independently until 
6^15 because the relative reaction rates are almost independent both on radius 
(except the {n, 7 ) reaction of and burnup. More exact investigations 

need consider that the relative reaction rates are not exactly constant as the 
burnup increasing. Such as from bu 2 to bug, Rg j decreases from 2.49 to 2.25, 
seen from Table 0] and O This is why the difference of isotope distributions 
between the MC and present calculations becomes larger as the increment 
of the burnup. One can use the present formula to obtain the concentration 
at a certain burnup level, such as bu^, re-simulate the reactions with the 
concentrations through MC, and refit the coefficients in the Eq. ([9]). 

The present calculation is limited to and ^^®Pu to examine 

its validity in a simple case. It is expected that it can be expanded to 
more complicated case, such as the consideration of other Pu isotopes and 
other nuclides which played important role in the reactor, such as the minor 
actinide and the poisons. In that case, a starting point is also needed to 
obtain the parameters in the Eq. ([9]). In addition, the present calculation 
constrains the Abu to compare with the MC calculations. In principle, the 
Abu can be rather small to obtain an almost continues changes of the reaction 
rates depending on the burnup and radius. In the case of small Abu, the 
assumption of constant reaction rates in Eq. ([7]) is more reasonable. 

For a more exact study, the present method can be used together with 
the MC calculation. Such as, starting from bui, the reaction rates and the 
concentrations can be calculated with rather small Abu up to bu 2 . The output 
is used for the next MC calculation. If the Abu is rather small, purely MC 
calculation is very time consuming from bui to bu 2 . But with the supplement 
of the present work, more exact and relatively quick results are expected. 

In practice, one can use the reaction rates in Eq. ([HD with the concentra- 
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tions at bui to calculate the AN{r) in Eq. ([7j) with small Abu, such as 0.01 
MWd/kgU. The next step is to use the new concentrations to calculate the 
reaction rates in Eq. (|9]) and then the next AN{r) through Eq. ([7]). After 
285 steps, the concentrations at bu 2 are obtained and can be used for the 
next MC calculation. It is actually the Euler method with very small steps, 
which is expected to be more exact. In Eq. ([7]), the reaction rates (actually 
the neutron flux and the concentrations) are assumed to be constant during 
Abu. But both the neutron flux and the concentrations change when the 
burnup changes. For detail, the concentration is in the reaction rate in the 
numerator of the Eq. ([7]) . The effect of the neutron flux in Eq. ([7]) is reflected 
by the integral of the reaction rates (the total hssion reaction rate) in the 
denominator. At the different burnup level, the total hssion rate and its 
distribution in each isotopes are different, which corresponds to the different 
neutron hux. If the calculation between bui and bu 2 are separated to many 
small steps, the neutron hux and concentrations are reconsidered for each 
steps. 

The present method is also helpful for solving the properties of the fuel 
rods at the diherent positions in an assembly or an reactor core. As discussed 
before, the relative reaction rates are almost constant in the fuel rod at the 
diherent position in the assembly. Only one factor is needed to describe the 
reaction rates in each fuel rod. One can simulate at an initial level through 
the MC method to obtain the factors for each fuel rod and investigate the 
evolution of the local burnup and isotope distributions for all fuel rods at 
certain average burnup level. 

4. Summary 

In conclusion, an analytical and simple formula is suggested to calculate 
the reaction rates as well as the burnup and the isotope distributions as the 
function of radius. The parameters of the formula are htted to the reac¬ 
tion rates of a given burnup level. Starting from the same burnup level, the 
present formula can give very nice description compared with a MC calcula¬ 
tion from the code TRIPOLI-4. The reaction rates depend on the concentra¬ 
tions and the neutron hux, which both vary along the radial direction and at 
the diherent burnup levels. The present work hnds almost constant relative 
reaction rates on these two degrees of freedom (except the {n, 7 ) reaction of 
on the degree of freedom of the radius), which provides a solid physical 
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explanation on the simple formula used to calculate the local burnup and the 
isotope distributions. 


Conflicts of Interest 

The authors declare that there is no conflict of interest regarding the 
publication of this manuscript. 

Acknowledgement 

The authors acknowledge to our collaborator CEA for the authorization 
of the TRIPOLI-4, and the useful discussions with Sicong Xiao, Min Zhang, 
Youjun Fang, and Nianwu Lan. This work has been supported by the Na¬ 
tional Natural Science Foundation of China under Grant No. 11305272, the 
Specialized Research Fund for the Doctoral Program of Higher Education un¬ 
der Grant No. 20130171120014, the Guangdong Natural Science Foundation 
under Grant No. 2014A030313217, and the Pearl River S&T Nova Program 
of Guangzhou under Grant No. 201506010060. 

References 

[1] K. Nogita, K. Une, J. Nucl. Mater. 226, 302 (1995). 

[2] C. B. Lee, Y. H. Jung, J. Nucl. Mater. 279, 207 (2000). 

[3] Y. H. Koo, et ai, J. Nucl. Mater. 295, 207 (2001). 

[4] J. Noirot, L. Desgranges, and J. Lamontagne, J. Nucl. Mater. 372, 318 
(2008). 

[5] M. Amaya, J. Nakamura, and T. Fuketa, J. Nucl. Mater. 392, 439 (2009). 

[6] V. Di Marcello, et al, Nucl. Eng. Des. 248, 149 (2012). 

[7] C. X. Yuan, et al, Nucl. Phys. A 883, 25 (2012); C. X. Yuan, et al, Phys 
Rev. C 85, 064324 (2012); C.X. Yuan, et al, ibid 89, 044327 (2014). 

[8] K. Lassmann, et al, J. Nucl. Mater. 208, 223 (1994). 

[9] K. Lassmann, C. T. Walker, J. van de Laar, J. Nucl. Mater. 255, 222 
(1998). 


18 



[10] A. Schubert, et al, J. Nucl. Mater. 376, 1 (2008). 

[11] C. B. Lee, et al, J. Nucl. Mater. 282, 196 (2000). 

[12] A. Soba, et al, J. Nucl. Mater. 433, 160 (2013). 

[13] A. Soba, et al, Ann. Nucl. Energy 70, 147 (2014). 

[14] M. Lemes, A. Soba, A. Denis, J. Nucl. Mater. 456, 174 (2015). 

[15] http://www.nndc.bnl.gov/exfor/endfOO.jsp 

[16] H. Matsumoto, M. Ouisloumen, T. Takeda, J. Nucl. Sci. Tec hn o). 42, 8 
(2005). 

[17] O. Petit, F. X. Hugot, Y. K. Lee, C. Jouanne, "TRIPOLI- 
4 version 4 user guide,” CEA-R- 6169, CEA (2008). 
http://www.nea.fr/abs/html/nea-1716.html (2008). 


19 




